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Peak Luminosities of Bursts from GRO J1744-28 measured with the 
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GRO J1744-28, discovered by BATSE in December 95, is the second neutron star system known to produce 
frequent accretion powered bursts. The system has been regularly monitored with the RXTE PCA since the peak 
of the first outburst in January 96 at which time the observed persistent and bursting count rates were ~ 25,000 
ct/sec and ~ 150,000 ct/sec, with corresponding instrumental deadtimes of ~ 10% and > 50%. We present a 
model which allows the reconstruction of the true incident count rate in the presence of enormous deadtime and 
use the model to derive a history of the peak luminosities and fluences of the bursts as a function of time. During 
the peak of the January 1996 and January 1997 outbursts, when the persistent emission was > 1 Crab, we infer 
peak luminosities of ~ 100 times the Eddington luminosity, and a ratio of persistent emission to integrated burst 
emission of ~ 34. 



1. GRO J1744-28 

GRO J1744-28, the "bursting pulsar" was di- 
covered by the Burst And Transient Source Ex- 
periment (BATSE) in December 1995 0] just 
prior to the launch of the Rossi X-ray Timing Ex- 
plorer (RXTE) on 30 December 1995. RXTE has 
performed a series of observations of this object 
from January 1996 through November 1997 dur- 
ing which time the source has varied from having 
a persistent flux in excess of 1 Crab with bursts 
which saturated the Proportional Counter Array 
(PCA) to being undetectable. The PCA observa- 
tions complement the nearly continuous BATSE 
record by virtue of the large collecting area and 
high time resolution of the PCA and High En- 
ergy X-ray Timing Experiment (HEXTE) exper- 
iments. These pointed experiments can study in- 
dividual bursts in much greater detail and follow 
the evolution of the outbursts to much fainter 
fluxes than BATSE. Studies of the individual 
bursts has, until now, been handicapped by an 
incomplete understanding of the deadtime pro- 
cesses in the PCA when observing extremely 



bright sources. The three main purposes of this 
contribution are to (a) present a model of dead- 
time processes in the PCA which is reliable for 
input fluxes at least up to 35 Crab ( ^ 0.5 x 10^ 
count/sec); (b) establish that the bursts from 
J1744-28 reached luminosities of ~ 100 times 
the Eddington luminosity; and (c) to provide ad- 
vice to users of RXTE who may observe sources 
brighter than 5 Crab during the rest of the RXTE 
mission. 

2. EVIDENCE OF SATURATION 

Figure 1 shows, for one burst observed on 1996 
February 4 and which peaked at 13:15:33 UTC, 
the observed good event rate as well as the Re- 
maining Counts rate from the Standard 1 data 
plotted with 0.125 second time bins. The Stan- 
dard 1 data (present for all PCA observations) 
contains 8 rates, read out every 0.125 seconds 
and internal calibration spectra. The eight rates 
consist of the good rate in each of the 5 detec- 
tors, the sum over 5 detectors of all events which 
trigger only the propane layer, all events which 
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Figure 1. Good and Coincidence rates observed 
from a 1996 February 4 burst. 



trigger the VLE flag (i.e. saturating events), and 
all other events commonly called the Remaining 
Counts rate. For most observations, the remain- 
ing Counts rate is dominated by particle induced 
background events; for bright sources there is a 
non trivial possibility that two (or more) cosmic 
X-rays will be detected in different parts of the de- 
tector and be recorded as a multiple anode event. 
The sudden rise in the Remaining Counts rate, 
coincident in time with the rise in the Good event 
rate, makes it clear that in the burst, the remain- 
ing counts rate is dominated by the coincidences 
of photons which arise in the burst. 

It is however obvious that while the good event 
rate reaches an apparent plateau around 1.5 x 10^ 
count/sec, the remaining count rate has a struc- 
ture with greater contrast. This observation sug- 
gests the central theme of our dcadtimc correc- 
tion: the rate of coincident events, which depends 
on the square of the incident rate, is a better mea- 
sure of the flux of very bright sources than the 
good rate. Information on the coincidence rate is 
always available on 0.125 second scales, and can 
be available with higher time resolution through 
the use of EDS modes which count and teleme- 
ter events which trigger exactly two anode chains. 
(For a complete list of modes see Q.) 

Examination of numerous bursts from January 



and February 1996 show similar characteristics: 
the good rate saturates near 1.5 x 10^ count/sec 
while the remaining rate shows much larger rela- 
tive variations. 

3. PCA AND EDS CHARACTERISTICS 

3.1. PCA 

Each Proportional Counter Unit (PCU) has 9 
independent signal chains. Seven signal chains, 
designated LI, Rl, L2, R2, L3, R3, and VP, 
signifying the "left" and "right" halves of the 
three Xenon layers and the propane veto layer 
have individual charge sensitive pre-amplifiers 
and shaping ampliflers and share a common ana- 
log to digital converter (ADC). The xenon veto 
anodes and the calibration flag set discriminators 
but are not pulse height analyzed. A good event 
is one which triggers exactly one signal chain; the 
pulse height from the ADC is unambiguously as- 
sociated with that signal chain. An event which 
triggers two or more chains produces two flags 
and a pulse height which could correspond to ei- 
ther chain (but in any case which contains only a 
fraction of the deposited energy) and is generally 
discarded as a non X-ray event. (An event which 
triggers the calibration flag and one other anode 
chain is an exception to this rule and is tagged as 
a calibration event.) The analog pulse shaping is 
a paralyzable process while the analog to digital 
conversion is a non-paralyzable process. 

3.2. EDS 

The EDS always runs two Standard modes in 
addition to the data modes selected by the ob- 
server. A key feature of each of these modes is 
that every PCA event is recorded exactly once. 
Standard 1 produces 8 rates read out every 0.125 
seconds while Standard 2 produces a spectrum 
for each signal chain and 29 distinct coincidence 
rates for each PCU once every 16 seconds. In spe- 
cial circumstances, the Standard 2 mode can be 
read out on shorter time scales. We present, in 
the next section, data obtained while observing 
bright sources with the Standard 2 mode being 
read out more frequently than usual. These data 
allow us to parameterize a model of PCA dead- 
time processes. 



4. DEADTIME MODEL 

Here we present data obtained while slowly 
scanning over Sco X-1 (with a 1 second readout 
interval for the Standard 2 data) and data ob- 
tained during a bright burst from J 1744-28 (with 
2 second readout). The details of the paralyzable 
deadtime process associated with the analog pulse 
shaping are pulse height dependent, so we expect 
some differences in parameters derived from Sco 
X-1 (a very soft source) and J1744-28 (a very hard 
source). 

Figure shows two rates observed during a scan 
over Sco X-1 (97Marl5). The smooth curve repre- 
sents the collimator efficiency for one of the PCU, 
and the histograms represent the counting rate 
observed in the LI layer (i.e. good events) and 
the rate of LI plus Rl coincidences (two photons 
observed within the coincidence window, one on 
each half of the first layer). The collimator effi- 
ciency is scaled to demonstrate that the LI rate is 
affected by a 25 — 30% deadtime near the peak of 
the collimator transmission. We believe that we 
can treat Sco X-1 as a constant source for these 
purposes; the presence of a nearly stationary 6 Hz 
QPO at the time of these observations indicates 
that Sco X-1 was in the normal branch with a 
momentarily stable mass accretion and luminos- 
ity. 

Figure || shows the same rate data, plotted 
against each other. The rate data have been cor- 
rected for the non paralyzable deadtime associ- 
ated with the analog to digital conversion (and 
which must account for all events in this detec- 
tor, not just the two rates that are plotted). Also 
shown is the best fit which models the Ll-I-Rl 
rate as a constant plus a term proportional to 
the square of the LI rate. The fit considers only 
points with a corrected LI rate < 10^ count/sec, 
although the figure shows the fit extrapolated to 
the highest observed rates. The constant repre- 
sents a background term, while the quadratic co- 
efficient represents the average time window dur- 
ing which a coincidence can be recorded. The 
average is over pulse height as we have ignored 
all energy information while creating this corre- 
lation. Identifying the LI and Rl rates as the 
same (due to their identical construction) , we ex- 
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Figure 2. Collimator transmission, LI rate, and 
Ll+Rl rate for one PCU during scan over Sco 
X-1 in normal branch 



pect that the quadratic term represents twice the 
coincidence window. The chance of an Ll+Rl 
event must be equal to the LI rate times the co- 
incidence window time the Rl rate, and we must 
multiply by two to account for Rl+Ll events, 
which are indistinguishable to the EDS. In this 
manner we identify the LI coincidence window 
6ti as 5.02/isec. 

We use similar fits to estimate two additional 
time windows: St2,Stp which are the coincidence 
windows associated with the second or third 
xenon layers and the propane layer, and a win- 
dow that describes the chance of getting an un- 
flagged event 6to. Unflagged events occur for 
events with certain separations on the same signal 
chains B . The process is very slightly more com- 
plicated when we fit Vp-|-Ll or Vp-|-Rl vs the LI 
rate. Here we assume that the window for having 
an LI event followed by a Vp event is known: the 
window for a Vp event followed by an LI event is 
then the quadratic coefRcient minus Sti. 

The entire process was repeated using the 2 sec- 
ond readouts and a number of bursts from J 1744- 
28 (96Jan27). Because the incident spectrum is 
quite different, the derived coincidence windows 
are different, particularly on the first layer. Table 
1 contains all of the derived windows. 



Ll+Rl = 5.3 + 1.033 • 10"= LI 
fit to LI < 10* 
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Figure 3. Ll+Rl rate fit as a quadratic function 
of LI rate. Both rates have been corrected for 
ADC deadtime. The fit is restricted to LI < 10'* 
but provides a reasonable description of the data 
at all observed rates. 
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Two additional pieces of information are 
needed to infer incident rates from observed co- 
incidence rates. First, the nonparalyzable dead- 
time associated with the analog to digital conver- 
sion is taken to be 9/zsec [Q. This value is inde- 
pendent of pulse height, and can be used for all 
sources. Second, we need the ratio of source re- 
lated events which interact in each layer of the 
detector, which is obviously dependent on the 
spectrum. For J 1744-28 we rely on the fact that 
the bursts apparently have the same spectrum as 
the persistent emission 13], and use the persis- 
tent emission to derive the ratios L1:L2:L3:VP 
to be 1000:150:75:180. (The rate on Rl is as- 
sumed equal to that on LI, and similarly for the 
other xenon layers). For Sco X-1, we measure 
the relative incident ratios from the scanning data 
obtained relatively far from the peak and obtain 
1000:67:21:440 (confirming our earlier statement 
that Sco X-1 has a much softer spectrum than 
J1744-28). 

4.1. Predicting the Coincidence Rate 

Our convention is that "incident rate" means 
the incident rate on the xenon layers. We can 
calculate the corresponding number of events ob- 
served in the propane layer. For a given incident 
rate, we calculate the incident rate on each signal 
chain Rj where the index j runs from 1 to 7 and 
corresponds to LI, Rl, L2, R2, L3, R3, and VP. 
The observed rate of 2 lower level discriminator 
events, R2LLD observed is given by 



R2LLD — AT^jY.i^jRjStjRi 



(1) 



where Stj is defined as 5ti for j = 1,2, as 6t2 
for j — 3,4,5,6, and as Stp for j — 7- A ac- 
counts for deadtime due to the non-paralyzable 
ADC process and is calculated as rm/(l -|- rintd)- 
Since these are predicted rates, we have complete 
knowledge about the assumed incident rate ri„. 
The observed rate of unflagged events, Rnofiag, 
and the observed rate of triply flagged events, 
RsLLD, are similarly calculated: 

Rnofiag — A'SjRjStoRj (2) 



RsLLD — ATijY^i^jTikyiijRjStjRiRk 



(3) 



In order to predict the Remaining Count rate as 
recorded in the Standard 1 data, all three terms 
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Figure 4. Good count rate, Coincident count rate, 
and two different estimates of the incident count 
rate. See text for details. The lower estimate of 
the corrected rate is our best estimate. 



must be added (for the bursts from J1744-28, the 
four fold coincidence rate shows a detectable rise 
associated with the burst, but it is not a signif- 
icant contributor to the remaining rate). In all 
cases, A must be evaluated considering all events 
(in other words, propane events always contribute 
to the ADC deadtime). 

Figure || shows the same data as figure ||, along 
with the corrected count rate. We corrected the 
count rate twice, once using the coefficients de- 
rived from J1744-28 and once using the coeffi- 
cients derived from Sco X-1. In both cases we 
assume the relative count rate distribution be- 
tween the layers appropriate for J1744-28. The 
higher correction, which peaks at about 4.6 x 10^ 
count/sec is associated with the coefficients for 
Sco X-1, while the curve that peaks near 4.2 x 10^ 
count/sec was derived using the J1744-28 co- 
efficients. That the two curves differ by only 
10% demonstrates the expected robustness of this 
technique if one or the other set of coefficients is 
used for a different bright source such as the Soft 
Gamma Repeater SGR 1806-20. It is easy to see 
why the Sco X-1 coefficients result in a higher in- 
ferred incident rate. Since we use the coincidence 
rate to infer the incident rate, since the count 
rate on the LI and Rl layers is the most impor- 



tant contributor, and since the coincidence win- 
dow derived from J1744-28 data is longer than 
that derived from the Sco X-1 data, we require 
a slightly lower incident rate to create the same 
double event rate. While 10% probably is a good 
estimate of the net uncertainty, we note that we 
have reconstructed this estimate in the presence 
of 70% deadtime by traditional definitions! Our 
estimate is similar to an independent estimate Iq] 
which derives relationships between the various 
coincidence rates and the good event rate without 
a detailed accounting for the different coincidence 
windows on different layers. While that model 
produces a similar result, and is therefore of equal 
phenomenological validity, they derive the value 
of the ADC deadtime, and consistently get a value 
which is larger than the measured 9/L/,sec. 

5. BURST HISTORY OF J1744-28 

Finally we apply the same deadtime correction 
discussed above to the data base of J1744-28 ob- 
servations. RXTE has observed J1744-28 regu- 
larly through out the first two years of the mis- 
sion (g g). Figure H shows the persistent flux, 
the peak ffiix in each burst, and the burst flu- 
ences. We have not included any bursts from 
May or June 1996; bursts at this time are ei- 
ther Type I X-ray bursts from another source 
within the field of view H or much fainter and 
more frequent bursts which are part of a different 
"rumbly" phenomenology B. The bright bursts 
occur only during the periods when the persis- 
tent flux is quite bright (above ~ 1000 count/sec). 
However, we have observed bursts over a range of 
persistent intensities that span one and a half or- 
ders of magnitude, and the striking feature of fig 
H is that the ratio of peak flux to persistent flux, 
and the ratio of burst fluence to persistent flux, 
remain relatively constant. 

Throughout the periods of bright bursting 
emission, the persistent emission can be ade- 
quately described by a hard power law modified 
by galactic absorption and a high energy expo- 
nential cutoff. (There is also evidence for an Iron 
line but this does not have a large effect on the de- 
rived flux). This simple model provides the same 
estimated flux as the continuous model described 
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Figure 5. Persistent flux (squares), peak burst 
flux (crosses), and burst fluence (circles) for burst 
observed during observations of J1744-28. Several 
bursts defy the general trend; these are largely 
Type I bursts from other sources in the field of 
view ([8]). 



by H. We find that the 2-10 keV fiux is 1.81 x 
10~^^ergsec^^ cm^^ per count/sec (counts/sec 
measured in all PCA channels); the 2-40 keV flux 
is 4.95 X 10~"'^^ergsec^^ cm^^ per count/sec; the 
2-100 keV flux is 5.40 x lO^^^ergsec"^ cm~^ per 
count/sec. Taking the highest estimated peak 
fluxes for bursts in the first outburst, with peak 
rates of 4.4 x 10^ count/sec, we derive a peak flux 
erg sec 



of 2.4 X 10"'^ 



^cm~^. The corresponding 
peak luminosity is 1.8 x lO^^d^ergsec"^ where ds 
is the distance to J1744-28 normalized to 8 kpc. 
This is about 100 times the Eddington luminosity 
given typical values for a neutron star mass and 
radius. Given the nearly constant ratio of burst 
fluencce to persistent flux, we conclude that the 
ratio a of persistent flux to time averaged burst 
flvix must be approximately constant (to the ex- 
tent that the average interval between bursts is 
constant). We previously reported that a ~ 34 
in January 1997 |1^] and take this as a represen- 
tative value for both major outbursts. 

None of the above conclusions could be reliably 
drawn from the good counting rate alone, but re- 
quire the extra information contained in the co- 



incidence rate. The most important conclsuion 
may therefore be our recommendation to future 
users of the RXTE satellite who may be observ- 
ing extremely bright sources: use modes which 
include some coincidence information at the high- 
est timescales of interest. On a practical level, ex- 
tremely bright sources may be thought of as those 
which produce > 15, 000 count/sec/PCU with as- 
sociated deadtimes of > 15%. An example of the 
value of this approach for observations of Sco X-1 
is given in [p|. 

Comparisons to the Crab nebula are often 
ambiguous due to the different spectral shapes. 
However, limiting the comparison to total count 
rate (13,000 count/sec for the Crab), we find that 
the bursts reached 34 Crab. Put another way, 
we can recast our subtitle (using "!" for "i"): A 
grand one man band with major lighting: stops at 
11 -^ Crah U 

two 
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